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SOLUTIONS OF BLUNT-BODY STAGNATION- REGION FLOWS 


WITH NONGRAY EMISSION AND ABSORPTION OF RADIATION BY 
A TIME -ASYMPTOTIC TECHNIQUE 

By Linwood B. Callis 
Langley Research Center 

SUMMARY 

A second- order time- asymptotic solution to radiation- coupled stagnation- region 
flows is presented. The solution is applied to the hypervelocity flow over blunt vehicles 
of inviscid, nonconducting, equilibrium air, emitting and absorbing nongray radiation. 
Velocities, nose radii, and altitudes covered by the analysis are sufficient to bracket 
reentry trajectories of current interest. Radiative heat-transfer rates for the range of 
interest and typical profiles of pressure, density, enthalpy, temperature, and velocity are 
shown. The nature of the time -asymptotic solution is discussed and it is shown to be a 
feasible means of achieving second-order accurate solutions to radiation- coupled shock- 
layer flows. 

Step-function models of the absorption coefficient are used in order to evaluate the 
divergence of the radiation flux vector. An analysis is carried out to determine what 
effect variations in the spectral complexity of the step model absorption coefficients used 
in the analysis will have on the thermodynamic and flow profiles of interest and on the 
nongray radiative heat- transfer rates. In this connection use is made of consistent model 
absorption coefficients having one to nine spectral steps with free-free, free-bound 
(including atomic line transitions), and molecular transitions taken into account. Rela- 
tively simple models of the absorption coefficient can be used with no significant loss of 
accuracy. An existing correlation for the cooling factor, the ratio of the radiation heat- 
transfer rate to the adiabatic radiation heat-transfer rate, is extended to larger veloc- 
ities than heretofore considered. 


INTRODUCTION 

During the past decade, problems involving radiation- coupled flow fields have 
received the attention of numerous investigators. It has been amply demonstrated (refs. 1 
to 4) that the inclusion of the emission and absorption of radiation in stagnation flow-field 
analyses may significantly alter both the flow-field structure and the heat-transfer rates 
typically expected from radiationless solutions. For blunted vehicles entering the earth's 



atmosphere, these effects manifest themselves at velocities in excess of 9 km/sec; this 
makes necessary the consideration of the coupling existing between the radiation trans- 
port equations and the fluid-dynamic equations governing the flow. The problem posed by 
such a consideration is extremely difficult since the solution is governed by a set of non- 
linear, partial, integro-differential equations. Furthermore, the integral terms involve 
complicated functions of wavelength and spatial coordinates; therefore, tedious numerical 
integration is necessary. 

The most general approaches to the problem to date are typified by the solutions of 
Howe and Viegas (ref. 1), Hoshizaki and Wilson (ref. 2), and Olstad (ref. 3). The solu- 
tions in references 1 and 2 are similar in that both consider the flow field to be totally 
viscous and conducting and both consider injection of foreign species as well as energy 
transport by diffusion. Howe and Viegas (ref. 1), assuming similarity, reduce their 
governing equations to a set of ordinary differential equations which are solved by stan- 
dard techniques. The absorption coefficient is considered to be gray. Hoshizaki and 
Wilson (ref. 2) maintain variations in the body tangent direction. The energy equation is 
solved by finite-difference techniques with the velocity and species concentration profiles 
determined from an integral solution to the species continuity equation and the momentum 
equation. The analysis is for a nongray gas. Olstad (ref. 3) uses a Poincare- Lighthill- 
Kuo singular perturbation technique to obtain stagnation region solutions for an inviscid, 
nonconducting gas with a nongray absorption coefficient. This solution involves expan- 
sions in terms of the radiation cooling parameter 46 a e and is valid for 45 a e « 1. 

The purpose of the present analysis is threefold. The first is to point out that time- 
asymptotic techniques are a practical means of solving nonadiabatic radiation- coupled 
flow fields. The techniques discussed herein have been used to generate radiation- 
coupled stagnation-line solutions, including emission and absorption of nongray radiation, 
and hopefully they will be applicable to the analysis of the full flow field. This particular 
approach was adopted as it seemed to offer a direct solution to the problem at hand with 
a minimum number of restrictive assumptions and computational difficulties. This 
advantage results in large part from the hyperbolic nature of the governing set of partial 
differential equations when cast in their unsteady form. Moretti and Abbett (ref. 5) have 
demonstrated that radiationless calculations can be carried on in the entire flow field 
with few of the usually attendant problems of convergence, stability, or sensitivity to ini- 
tial guesses. The time- asymptotic approach used herein is essentially a second-order 
accurate Taylor series expansion in time of the quantities of interest. The solution is 
advanced in time until an asymptotic situation is reached representing the steady- state 
result. 

Shock waves, with the aid of unsteady characteristics, are treated as discontinuities 
providing more accurate profiles in the vicinity of the shock front. In addition, use of a 
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discontinuous shock allows a coarser mesh to be used than is possible with the "artificial 
viscosity" technique (refs. 6 and 7); thereby, there is a considerable saving in computing 
time and storage. The second-order solutions presented in this report have been 
achieved in remarkably short computing times, typically 20 seconds/case or less on a 
Control Data series 6000 computer system. The present stagnation solution is thus 
offered as a demonstration of the feasibility of the time-asymptotic approach to the solu- 
tion of general radiation- coupled flow fields as well as a valid solution in its own right. 

The second purpose of the analysis is to determine what effect variations in the 
spectral complexity of the absorption coefficients used in the analysis will have on the 
thermodynamic and flow profiles of interest and on the nongray radiative heat-transfer 
rates. In this connection, use is made of data tabulated by Olstad (ref. 8) in formulating 
step models of the absorption coefficient with one to nine spectral steps. The radiative 
processes taken into account in the formulation of these step models are line transitions, 
N + - and 0 + - electron recombination, free-free electron interactions, and molecular -band 
system radiation. Solutions determined with model absorption coefficients constructed 
from the same detailed spectral information are compared and discussed. 

Finally, radiation heat-transfer rates to a sphere are presented for velocities from 
9 to 18 kilometers/second, nose radii from 0.1 to 10 meters, and altitudes from 49 to 
73 kilometers. A comparison is made with the cooling factor correlation presented in 
reference 8, and the range of this correlation is extended, 

SYMBOLS 

A parameter defined by equations (A18) and (A23) 

B,C,D,E,Pj,P 3 parameters defined by equations (A18) 

B„ nondimensional blackbody function 

Cj,C 2 ,C 3 parameter defined by equations (B4), (B5), and (B6), respectively 

C °° — Yt 

E n exponential integral function of order n, E n (Y) = \ e t -n dt 

J 1 

Fr nondimensional divergence of the radiation flux vector 

F r cooling factor, 

C ^R,w,a 

f parameter defined by equation (A24) 
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parameter representing p, p, v, or f 
nondimensional enthalpy 
nondimensional specific radiation intensity 
parameter defined by equation (A10) 
Boltzmann constant, joules/°K 
nondimensional pressure 
nondimensional radiation heat flux 
radius defined by equation (A5) 
nose radius, meters 

nondimensional incremental radiation path 
nondimensional time except where indicated 


T* parameter defined by equation (C6) 

U nondimensional velocity vector 

u,v nondimensional velocity components 

Voo free-stream velocity, meters/second 

x’,y',z' nondimensional Cartesian coordinates 

x,y nondimensional body- oriented coordinates 

X,Y nondimensional floating coordinates 

Z altitude, kilometers 

a^,a parameters describing picket-fence model 
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Subscripts: 

a 

o 

w 

A 


angle between shock and body except where indicated, radians 

nondimensional standoff distance 

nondimensional time increment 

parameter defined by equation (A3) 

parameter defined by equation (C7) 

body surface inclination, radians 

wavelength, angstroms except where indicated 

nondimensional radiation frequency 

density, kilograms/ meter3 

absorption coefficient, meter" 1 

Planck mean absorption coefficient, meter" ^ 

Stefan- Boltzmann constant, watts/ (meter-OK^) 
optical thickness 
azimuthal angle, radians 
solid angle, steradians 

adiabatic wall 

standard atmospheric conditions 
wall conditions 

stationary shock conditions on center line 
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local unsteady shock conditions at time T 


v quantity evaluated at frequency v 

°° free -stream conditions 

Symbols with bars represent dimensional quantities, those without bars indicate 
nondimensional quantities. 


ANALYSIS 


Nondimensionalization 

The equations appearing in this report are nondimensional. All distances are non- 
dimensionalized by the vehicle nose radius Rjq-, all time quantities by the ratio Rjq/ V oo> 

densities by p^, pressures by p^sf 0 a , energy flux by ip^V , enthalpies by Voo , the 

divergence of the radiation flux vector by a^pkp^T^^ where is a temperature, 

velocities by Voo, and absorption coefficients by (pk p ^. 

Basic Assumptions and Flow Model 

The basic assumptions under which the analysis is simplified and carried out are 

that 

(1) The flow is axisymmetric, inviscid, and nonconducting 

(2) The flow is in chemical equilibrium 

(3) The shockwave and free stream are transparent to shock- layer radiation 

(4) The body surface with regard to radiation is cold, nonreflecting, and black 

(5) The tangent slab approximation is applicable 

Basically these assumptions are made so that the applicability of the time- dependent tech- 
niques to radiation- coupled flows can be clearly determined free of other complicating 
phenomena. They are discussed at various points in the analysis. Figure 1 presents a 
schematic illustration of the flow geometry in the stagnation region of spherically blunted 
vehicle. 


Procedure 

The time- asymptotic approach in this report uses the following second- order accu- 
rate Taylor series expansion to advance an assumed solution in time until it becomes 
asymptotic to the desired steady- state solution: 



( 1 ) 


g(Y,T + AT) = g(Y,T) + g T (Y,T)AT + | g TT (Y,T)AT 2 + O^T 3 ) 


The expansion furnishes a new value of the parameter g, which may represent p, p, 
v, or f, at the coordinates Y and T + AT. This procedure is carried out at each Y 
coordinate of interest until all the information required for the continuation of the solu- 
tion in time is available at T + AT. Information on the shock and the body at T + AT 
is generated with the aid of the method of unsteady characteristics, discussed in the next 
section. Once the information at the next time plane is completely determined, the pro- 
cess is repeated until the variation between solutions at successive time planes is suffi- 
ciently small. This is taken as the steady-state result. 


The first-order time derivatives are determined from the unsteady conservation 
equations developed in appendix A (eqs. (A20), (A21), (A22), (A23), and (A25)). These 
equations are given here as follows: 


. . 2pf py Y 2Kpv 
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The right-hand side of equations (A20), (A21), (A22), (A23), and (A25) (hereafter referred 
to as set 1) are evaluated by using data on the most recent time plane. All derivatives, 
with the exception of p-^, are determined by centered finite difference representation. 

At this point, it is necessary to determine values for p-^. This is accomplished by 
assuming that at the center line the shock wave is concentric with the body so that 

p XX = " 2 ( 1 ‘ pg) ^ 

Equation (2), derived from steady-flow shock relations, is admittedly inconsistent with 
the unsteady nature of the analysis. The consequences of this inconsistency are not 
serious, however, since the momentum equation is only weakly coupled to the energy 
equation and since a steady- state solution is ultimately sought. Hence, use of equation (2) 
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involves some error in the intermediate development of the solution but little in the final 
results. 


Second- order time derivatives in equation (1) are obtained by differentiating set 1 
with respect to T. The mixed space-time derivatives that result are evaluated by dif- 
ferentiating set 1 with respect to Y and using centered finite differences where appli- 
cable. Time derivatives of quantities such as P 3 , P^, Fr, and A are cast as back- 

ward time differences by using the most recent time data and data from the previous time 
plane. The result of these manipulations is a set of equations for the required first and 
second time derivatives of p, p, f, and v. These derivatives used in equation (1) per- 
mit the solution to advance in time. 

Analysis of Shock and Body Points 

A quasi- one- dimensional unsteady characteristics technique, similar to that used 
by Moretti and Abbett (ref . 5) was applied to the analysis of the shock wave in order to 
retain the shock as a discrete surface. Artificial viscosity techniques have the effect of 
smearing the shock until it becomes not a discrete surface but a region having large flow 
property gradients. The distinction between the two approaches can, for radiation- 
coupled flows, be an important one since the enthalpy level can be significantly reduced in 
a region near the shock front. Consequently, poor shock definition can result in a poorly 
defined enthalpy profile. A schematic of the flow field (fig. 2) near the shock serves to 
illustrate how the characteristics method is applied. The shock standoff distance 5 
varies with T. At time T all quantities are known across the shock layer and it is 
desired to know 6(T + AT), 6 -j>(T + AT), and values of the flow properties at point A. 

As an initial guess, it is assumed that 6 T (T + AT) = 6 T (T) and that 
6(T + AT) = 6(T) + 6rp(T)AT; thus, point A is located and, with the aid of the Rankine- 
Hugoniot relations, conditions are determined at A in the shock layer. A right running 
characteristic (see appendix B), indicated in figure 2, is extended in the negative time 
direction until B is located. The location of B is stabilized by an iterative process aver- 
aging the slope of the right running characteristic between A and B. Once point B is 
stabilized, the compatibility relation (eq. (Bll)) is integrated from point B to point A with 
averaged coefficients and a value of v at point A determined. If the computed v does 
not match the value obtained from the Rankine - Hugoniot equations a new value of 
5'j'(T + AT), consistent with v and A as determined from the compatibility equations, is 
calculated and the entire procedure cycled to convergence. This procedure yields condi- 
tions at point A and is carried out prior to the advancement of interior flow properties 
to T + AT. 

At the body, the procedure is similar. The slope of the left running characteristic 
(fig. 2) is assumed at A’ and extended in negative time until point B’ is located. Point B’ 
is stabilized by iteration and the proper compatibility relation (eq. (Bll)) integrated from 
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B' to A' to determine the pressure at A'. The unsteady energy equation is then integrated 
from C* to A' to determine the enthalpy at A'. Values of h and p permit a new calcu- 
lation of the characteristic slope at A' and the procedure is cycled to convergence. 

Computation of the Divergence of the Radiation Flux Vector 

The calculation of the divergence of the radiation flux vector F R is carried out 
under the assumptions and by the means indicated in appendix C. It has been determined 
that F R need not be calculated for each time cycle but may be evaluated, held fixed for 
several cycles, and then recomputed. This saves considerable computational time and 
has a negligible influence on the final result. The transient portion of the solution obvi- 
ously will differ from that in which F R is determined anew each cycle. 

Step Model Absorption Coefficient 

The difficulty of and the time consumed in the calculation of the divergence of the 
radiation flux vector F R is strongly dependent on the spectral complexity of the absorp- 
tion coefficient used. If a true replica of the absorption coefficient including line, free- 
bound, and free-free transitions were used, the computation time for F R would be such 
that general flow-field calculations about blunted bodies would be impractical. As a 
result, various simplifications to the detailed absorption coefficient have been offered. 

The most obvious approximation and, also, the crudest is the gray approximation in which 
the absorption coefficient is allowed to vary with temperature and density but not with 
wavelength. The gray coefficient is usually the Planck mean of the true absorption coef- 
ficient. This approximation is unacceptable since it eliminates spectral detail which may 
be important when the shock layer is not optically thin in all spectral regions. A second 
approximate way to represent the absorption coefficient is to use two or more spectral 
steps to model the true coefficient. 

The height or value of pk for each of the steps varies with temperature and den- 
sity but not with wavelength over the spectral extent of the step. Clearly, with the addi- 
tion of a large number of steps, the coefficient may be approximated as closely as 
desired. 

The approach adopted herein is that, for engineering purposes, an adequate repre- 
sentation of the absorption coefficient can be realized with a step model. The two-step, 
three-step, and nine-step models of pk used in this report are discussed in appendix D 
and summarized in tables I and II. For comparison purposes a gray model is included. 

A schematic of these models is shown in figure 3. 
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Equilibrium Air Properties 

Expressions used in this report relating pressure and sound speed to p and h 
for equilibrium air are presented by Cohen (refs. 9 and 10) or have been derived from 
data presented in these references. For the temperature, correlations due to Olstad 
have been used. More exact thermodynamic data, either tabulated or calculated, were 
not used in an effort to save computational time and in the interest of program simplicity. 
With regard to the general accuracy of the overall solution, it is felt that uncertainties in 
the absorption coefficient information generally available outweigh the slight errors due 
to the use of thermodynamic correlations. 

Stability 

The stability requirement applied to the present analysis is the CFL (Courant- 
Friedricks- Lewy) criterion used by Moretti and Abbett (ref. 5). The point of view was 
adopted that no consideration (other than optimization of the time increment) would be 
given to stability analyses unless instabilities developed in the solution. Since no such 
instabilities were encountered during the development and utilization of the present solu- 
tion, the inclusion of a stability analysis is omitted. The application of the CFL criterion 
has proved satisfactory. 

The uniqueness of the present solutions was determined by the execution of a num- 
ber of different cases with significantly different initial conditions. In each case, varia- 
tions of the initial conditions caused the final result to be changed by less than 1/2 per- 
cent. This observation and comparisons with other solutions led to the conclusion that 
no uniqueness problems existed. 


RESULTS AND DISCUSSION 


The analysis described has been successfully applied to flow situations with 
9 i = 20 km/sec, 0.01 s R^j s 10 m, 36 i Z g 73 km. Typical results for p, p, v, 
T, and h are shown in figure 4. The variations in p and h are as expected for a 
radiation- coupled flow field. Figure 5 depicts the variation in the enthalpy profile for 
three values of V^. The enthalpy level is reduced as Vco increases because of the 
increased energy lost by radiation from the shock layer. A similar variation in the 
enthalpy level is apparent as the nose radius is increased. (See fig. 6.) Variation 

of the nondimensional standoff distance with both Rjj and Voo is shown in figure 7. At 
the smaller values of Rjj, 6 approaches its adiabatic value, whereas at larger values 
of Rjj, 6 is significantly reduced as a result of the higher density levels in the shock 
layer. 
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The results shown in figures 4 to 7 demonstrate that the present time- asymptotic 
technique can be used successfully in the analysis of radiation- coupled flows. Results 
presented in this report cover the widest range of conditions yet published for a single 
solution. Validation of these results is accomplished by comparison with the work of 
reference 8, which uses an absorption coefficient identical to that shown in table I, and 
by comparison with the work of several other authors using different absorption coeffi- 
cients. These comparisons are included in the discussion of figures 10 and 11, respec- 
tively. Typical computing times for the solutions are approximately 20 seconds/case 
or less (Control Data series 6000 computer system) making the technique attractive from 
the standpoint of practical usage. Computing times are relatively independent of Voo, 

Rn, and Z. 

Figure 8 shows, for a typical velocity and ambient density and a range of nose radii, 
the contribution to the total radiative heat flux by the various spectral steps of the nine- 
step model. For this ease, the equilibrium mole fractions of selected air species are 
given at the shock and body in table in (1% = 1 m )- Several observations can be made 
from figure 8, the first being that the largest single contribution to the total flux for all 
Rjq- comes from step 9, which has the smallest optical thickness. For R N £ 1.5 meters 
this contribution is larger than the sum of the remaining contributions, an indication of 
the effect of reabsorption on the contribution from the more optically thick steps. A 
second point of interest is the relatively modest contribution made by the line centers, 
steps 4, 6, and 8. Steps 4 and 6 contribute such a minor amount that they may be com- 
bined with other steps with very little error. Step 8, however, can be as much as 6 to 
10 percent of the total at small R^-. Its identity should be preserved. Steps 1, 2, 3, 
and 5 all contribute heavily at small Rjj. These steps, together with 4 and 6, are com- 
bined according to a Planck average to form a single step existing spectrally from 400 
to 1130 A. Two other steps may be formed by the combination of 7 and 9 which behave 
similarly and by the remaining step 8. The result is the three-step model indicated in 
table II. The two-step model of table n is formed with the inclusion of step 8 with 
steps 7 and 9. 

It should be noted here that the values of the absorption coefficient for steps 4, 6, 
and 8 are based on Griem's f-numbers (ref. 11, the best available when the model was 
formulated) which, according to more recent information, may be low by an order of 
magnitude. If updated values of the f-numbers were used, conclusions concerning the 
effect of the line centers on the heat flux may be modified. However, the same type of 
argument could be applied. 

Figures 9 and 10, for a typical flow condition, show the effect, on both the heat- 
transfer rate and the enthalpy profiles, of using the nine- step, three- step, two-step, and 
gray absorption coefficients. H the nine-step model is accepted as being most nearly 
representative of the exact absorption coefficient, it is clear that the gray model is 
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unacceptable except for qualitative analysis, the enthalpy profile being in error by as 
much as 30 percent. The three- step model has a maximum enthalpy error of approxi- 
mately 3 percent and the two-step model of 15 percent. This is an indication that the 
contribution of the line centers, given by step 8, should be maintained as a separate step. 

A large nose radius (in this case, 8 meters) emphasizes the defects of the two- and three- 
step models. In figure 10, stagnation-point heating rates are compared with those of ref- 
erence 8 over a range of and two values of Voo. It should be pointed out that the 

solutions of reference 8 use the nine- step absorption coefficient (table I), whereas the 
present solutions are for both the nine- step and the three-step models (table II). The 
solutions of Olstad (ref. 8) are terminated at values of Rn and Voo at which the 
approximations of the small, perturbation solution begin to be suspect. It is felt that some 
of the divergence between the solutions at the larger nose radii can be attributed to the 
onset of the breakdown of the small perturbation solution.. Over most of the range of 
Rjg, however, the agreement is excellent, slight differences in the level of the two nine- 
step models being attributed to the use of different thermodynamic correlations. The 
close agreement between solutions using the three- and nine- step models is interesting 
in light of the variation of the complexity of the step models used. 

In figure 11, the present stagnation-point heating rates are compared (for a range 
of Voo and P^/Po) w ^h those of Hoshizaki and Wilson (ref. 2), Wilson and Hoshizaki 
(ref. 12), Chin (ref. 13), and Page, Compton, Borucki, Cliffone, and Cooper (ref. 14). All 
these references, with the exception of reference 2, include the effects of line radiation. 
Figure 11, though not the best way to present a comparison of solutions, is convenient 
since most of the data included were readily available in this form. (See ref. 14.) As 
seen from figure 11, the data of Wilson and Hoshizaki (ref. 12) with p^/p 0 = 0.3 and 0.4 
compare very well with the present data for 0.08 < P^/P 0 = 0.2, the radiation heat- 
transfer rate being weakly proportional to P^/P 0 - The data of Chin and of Page and 
coworkers (refs. 13 and 14, respectively) appear to be slightly low over the velocity 
range shown. The data of Hoshizaki and Wilson (ref . 2) which take into account continuum 
radiation only are 41 percent below the present result emphasizing the importance of line 
transitions. In general, however, the comparisons shown in figure 11 are satisfactory. 

Figure 12 presents, for a spherically tipped body, a range of stagnation-point radia- 
tion heat-transfer solutions calculated by using the three- step absorption-coefficient 
model. Solutions are given for a range of R N and values of V M and Z sufficient to 
bracket most reentry conditions of interest. 

In reference 8, Olstad presents a correlation of the cooling factor F c , the ratio of 
the stagnation- point radiation heat-transfer rate q„ to the adiabatic radiant heat- 
transfer rate q R w a , both heating rates including the effects of absorption. The cooling 
factor is correlated as a function of q R w a for various altitudes and nose radii and for 
a fixed velocity. 
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Figure 13 presents these correlations for F c as they appear in reference 8. The 
present solution has been used to extend the velocity range to 18 280 m/sec and to extend 
the correlations to larger values of q*, For values of q D . . such that both solu- 

XV J W ^ d. llj VV y Or 

tions are equally applicable, excellent agreement is obtained. At the larger values of 
qR a , there is a tendency for the solutions to diverge, the present solutions having 
slightly higher values of F c . This is attributed to the aforementioned onset of the break- 
down of the small perturbation solution of reference 8. For further comparison, data 
from references 2 and 15 are indicated in figure 13(c), all falling well within the correla- 
tion limits. The data of Hoshizaki and Wilson (ref. 2) include no effect of line radiation, 
whereas the data of Rigdon, Dirling, and Thomas (ref. 15) use a modified hydrogenic 
approximation. In terms of absolute heating rates, agreement between the results of 
these references and the present results could differ by as much as 40 percent. (See 
fig. 11.) However, the good agreement of the data of references 2, 8, and 15 and the pres- 
ent data when correlated as in figure 13 is taken as a substantial indication of the validity 
of the correlation and its effectiveness in taking out the effect of variations in the absorp- 
tion coefficient. 

In light of these observations, it appears reasonable for an investigator to devote 
considerable effort to an accurate determination of q 0 „ (essentially a constant tern- 
perature calculation) with the best available absorption coefficient. The quantity qR 
may be translated with a high degree of accuracy (5 percent) to the desired heating rate 
qR . This procedure would eliminate many of the difficulties and much of the time asso- 
ciated with the fully coupled calculation. The determination of qR w could be further 
simplified by the compilation of tables of q„ for a range of conditions. This would 

itjWj ci 

be useful for design purposes. 

CONCLUDING REMARKS 

It is felt that the present analysis successfully demonstrates the applicability of 
time- asymptotic techniques to the solution of inviscid, nonconducting flows including the 
effects of emission and absorption of nongray radiation. The advantages of the method 
are that it is accurate to second order, that few approximations are necessary, that, little 
computational time is required, that solutions are insensitive to initial input, and that no 
stability problems are encountered if the CFL stability criterion is satisfied. With 
regard to full flow-field calculations, solutions are direct and may be easily extended 
through the sonic line with few of the usually attendant problems of other methods. 

The effect of the variation of the spectral complexity of the model absorption coef- 
ficients on the final solution profiles and heat-transfer rates has been investigated. It 
has been found that the three- step model can be used to adequately represent the nine- 
step model over a range of velocities and nose radii. The question raised at this point 
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is whether or not a simplified model of the absorption coefficient (such as the nine-step 
model or one more sophisticated) can be used in practical calculations and yield suffi- 
ciently accurate results (less than 15-percent error) when compared with calculations 
using detailed spectral absorption coefficients. Obviously the only rigorous way to answer 
the question is to calculate both divergence of the radiation flux vector and the energy flux 
throughout a nonisothermal slab using detailed spectral information. These calculations 
should then be compared to similar calculations using modeled absorption coefficients 
derived from the detailed spectral information. To the author's knowledge, this has not 
yet been done. It is felt, however, that the present analysis strongly indicates that the 
detailed absorption can be simply modeled and these models used to achieve results with 
good accuracy and the retention of all significant trends. If full flow-field calculations 
are to be carried out, such work would be greatly facilitated by the use of modeled 
absorption coefficients. 

A wide range of stagnation- point radiation heat- transfer solutions is presented 
which allows a direct determination of the heating rate for the nose radii, velocities, and 
altitudes shown. In addition, extensions of previously published cooling factor correla- 
tions have been made and are presented. These extensions allow the determination of 
the radiation heating rate for arbitrary altitudes, nose radii, and (within certain limits) 
absorption coefficients. Required for this determination is the velocity of the vehicle 
and a calculation of adiabatic wall heat flux by using the chosen absorption coefficient. 

Use of these correlations greatly simplifies the determination of the nongray radiation 
heat transfer for fully coupled flows and should prove valuable for design work. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., September 24, 1968, 

129-01-03-09-23. 
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APPENDIX A 


COORDINATE SYSTEM AND BASIC EQUATIONS 


Coordinate System 

The coordinate system used in the analysis is a "floating" system (sketch Al) 
obtained from a transformation of body-oriented coordinates (t,x,y) for axisymmetric 
flow (sketch A2). The procedure will be to transform the equations of motion from the 
Cartesian system (t',x',y’,z’) to the body-oriented system, then to recast these relations 
in the "floating" coordinates (T,X,Y). 



Sketch Al 
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Basic Equations 

The following unsteady inviscid, nonconducting equations, including radiation effects, 
are given as follows: 

Continuity 

|L + v.(pu) = 0 (Al) 

Momentum 


Energy 


— + TJ • VU - - I Vp 

at' p y 



6 F R 1 3p 
P p at’ 




(A2) 


(A3) 


The relations governing the transformation from the Cartesian coordinates (t',x',y',z') 
to the body-oriented system (t,x,y) are as follows: 


where 


t* = t 


x' = cos 9 dx - y sin 6 


y ' = r sin cj> 
z' = r cos 4 > 



sin 9 dx + y cos 9 


(A4) 


(A5) 


If the proper transformation equations are applied and, for axisymmetric flow, the 
azimuthal components ignored, equations (Al) to (A3) become, in the body- oriented system, 

Continuity 

f? + X 4 (pu) + ? sin 6 + + pv (j + £2 F^) = 0 (A6) 

x-momentum 


8u u 3u 8u Kuv _ 1 9p 

at X 8x + 8y + X “ pX 3x 


(A7) 
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y- momentum 

£Z. + H.^Z + v — - Ku ^ 1 9 P 

at A 9x 9y A P 9y 

Energy 

9h U 8h 8h _ 1 9p u 9p V 8p eF R 

9t A 9x 9y P 9t pA 9x P 9y " P 


(A8) 


(A9) 


where 


A = 1 + Kyi 



(AlO) 


It would be convenient, if during the unsteady calculation the shock itself became a 
coordinate line as well as the body. This is accomplished through the transformation to 
the "floating" coordinate system shown as follows: 



Y = l 


(All) 


which gives for y = 0, Y = 0 and for y = 6, Y = 1. The differential relationships 
between the two coordinate systems are given below 


9Q = 9Q_ Y 95 9( ) 
9t ST “ 6 9T 9Y 


9X 



9y 6 9Y 


(A 12) 


With the state equation, 


h = h(p,p) 


(A13) 


and switching to p as the dependent variable in the energy equation, equations (A6) to 
(AlO) under the transformation (All) become 
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Continuity 


X- momentum 


Y- momentum 


Energy 


p-p = 


= '( Ap Y + I p x- Bu y + x u x + S v y + c ) 


U T - ^Au y + “ u X ■ ^ p Y + Jx p x + d 


= '( Av Y + I t X + ^6 P Y- E ) 


Ap Y + ! p X + P 3(f “X - Bu Y + § V Y + C ) - 


eF 


R 


P 1 P 


(A 14) 


(A 15) 


(A16) 


(A17) 


Subscript notation has been used to denote partial differentiation. The parameters in the 
equations are defined by the following relations: 

A - - - _ HX. tan 3 — ^ 

A - g g 6 tan P x 


B = 


_ a 6 pY 


X 5 


tan /3 


C = ^ sin 0 + P v 


cos (A 


D = 


E = 


Kuv 


Ku^ 


P.3 


= 9h/8p 


1 

P 


ah 

ip 


P = — - 

1 P ap 


(A18) 


Equations (A14) to (A18) are valid for inviscid, nonconducting, equilibrium flows 
over smooth blunt bodies including emission and absorption of nongray radiation. 

For this analysis these equations must be reduced so that they are applicable to 
stagnation- line flows. Symmetry dictates the following conditions at the stagnation line: 


18 



APPENDIX A 


U = Uy = 0 

Px = v x = h x = P x = 0 


(A19) 


Indeterminant terms such as pu ^ in - and a re evaluated with the aid of 

L' Hospital’s rule. Equations (A14) to (A18) now reduce to the following set: 


Continuity 


Y- momentum 


P T - -(a p y + ^ u x + -^ + 


V T 


Vm = ~ AV„ + 


'Y + pS y Y 


Energy 


P T = 


Ap + P,(i£ u , | eF R 

AP Y + P 3\x U X + 6 + X ) P lP 


where under the stagnation line restrictions (eqs. (A19)), 

.. Y6 m 


A = 


5 


(A20) 


(A21) 


(A22) 


(A23) 


Equation (A15), the X- momentum equation, becomes degenerate at the stagnation line and 
seemingly offers no useful information. However, if equation (A7) is differentiated with 
respect to x, restricted to the stagnation line, and expressed in terms of the defined 
quantity f, where 


f = u x = u x 


(A24) 


the following differential equation for f is obtained: 

f /*f f 2 , Kfv , P XX\ / A9t .\ 

f T = - \Ai Y + - + — + —j (A25) 

after transformation to the (T,X,Y) system. This equation is necessary because of the 
appearance in equations (A20) and (A22) of ux- 

The governing differential equations are now (A20), (A21), (A22), and (A25). 
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DEVELOPMENT OF UNSTEADY CHARACTERISTICS EQUATIONS 
IN FLOATING COORDINATE SYSTEM 

A quasi-one-dimensional unsteady characteristics solution, similar to that reported 
in reference 5, is used for the shock and body points. The compatibility equations and 
characteristic directions are determined from the general equations cast in the following 
form. The asterisks in equations (Bl) and (B4) indicate that it is the natural logarithm 
of p that is differentiated.. 

P T * + A p Y * + | v Y = Cj (Bl) 

v T + Av y + i p Y - C 2 (B2) 

Pop 

p T + Ap Y + — g- v y = C 3 (B3) 

where 

C l = '(l p X* + X'| U Y + f) 

c 2 = -(f v x - E ) (B5) 

c 3 - -[! Px - 1 + c ) - - P 3 Bu r| 0M> 

At this point, it should be noted that, in the vicinity of the shock wave, the pertinent quan- 
tities are the normal component of the velocity and the Y-coordinate. The tangential 
velocity is unchanged across the shock and for smooth bodies the tangential derivatives 
of the flow properties are, for the most part, small in comparison with other terms. 
These observations coupled with the fact that for hypervelocity flows, the Y coordinate 
in the "floating" system is nearly normal to the shock (it is normal on the stagnation line 
for axisymmetric flow) suggest treating the terms represented by C^, C 2 , and C 3 as 
driving terms and employing a one-dimensional characteristics treatment of the shock 
point. The same basic argument can be made with regard to body points and it is for 
these reasons that equations (Bl) to (B3) are cast in their present form. A variation of 
this basic technique was discussed by Moretti in reference 16 and more detailed informa- 
tion can be found in that publication. 

To determine the characteristic directions for equations (Bl) to (B3) and the fol- 
lowing differential relations: 


20 



APPENDIX B 


^ dT + ^dY = dp* 

(B7) 

||dT + ||dY = dv 

(B 8 ) 


(B9) 


the associated coefficient determinant is set equal to zero and expanded. The resulting 
cubic in dY/dT is solved and the characteristic directions are found to be 


dY 

dT 



(BIO) 


The compatibility equations are determined by replacing a column of the coefficient 
determinant by the column formed by the right-hand side of equations (Bl) to (B3) and 
(B7) to (B9). Expanding the resulting determinant and substituting equation (BIO) yields 
the following equation: 


dv _±_ dp = c 
dT dT pj 


(Bll) 


Equations (Bl) to (Bll) are valid for general flow fields. For axisymmetric 
stagnation- line flow equations (BIO) and (Bll) maintain the same general form; however, 
A is given by equation (A23) rather than equation (A18), Cg = 0, and C 3 becomes 


C3 = 


2pP 3 

X 


(f + Kv) 


eF. 


R 


PlPj 


(B12) 
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EVALUATION OF THE DIVERGENCE OF THE RADIATION FLUX VECTOR 

It is desired to derive and determine a means of evaluating the expression for the 
divergence of the radiation flux vector at points across the shock layer. The shock layer 
at the stagnation line is modeled as an infinite gas slab of thickness 6 with one boundary 
transparent to incident radiation and nonemitting. The other boundary, representing the 
body surface, is to be a cold, nonemitting, reflecting surface. Properties are assumed 
to vary only with the normal coordinate, and local thermodynamic equilibrium is assumed. 

The nondimensional divergence of the flux vector is represented by the integral 

r 00 r 8l i> 

F r = \ A — — do? dv (Cl) 

R ^4tt 9s 

where l v is the specific intensity of the radiation, and ds an incremental length along 
the path a pencil of radiation travels. The radiative transfer equation is substituted and 
after integrating there results 

F r = 4(pk p )T 4 - (pk^Iy dw dv (C2) 

If Ij,, as determined from a solution of the radiative transfer equation, is substituted and 
the proper boundary conditions applied, equation (C2) becomes 

F r = 4^pkp^T 4 2tt ^ (pkj^ ^ - t^j^dt^, - 2 r w E 2 ^T^^ B t ,E 2 (t I ^dt {/ dv 

(C3) 

where r w is the constant coefficient of reflection of the wall and t v a dummy variable 
of integration. 

In the present report r w = 0, and the absorption coefficient, discussed in appen- 
dix D, is approximated by a step function model with one to nine spectral steps. The 
value of the absorption coefficient is invariant within spectral range of each step but is 
allowed to change with temperature and density. Under these assumptions equation (C3) 
may be written as 
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where N is the number of spectral steps considered. Since pkj is constant over 

it may be removed from under the integral and the integrations over v and t, the 
optical thickness, transposed. Equation (C4) then becomes 

N t 

T-I A / _l_ \ m4 o \ ' / *1. \ r* * ^4 x3 / f „ 4.. l\jl. t r* c\ 


Er ~ 4(pk p )T 4 - 2 £ (pki) J o ’V^E^n - t^dt* 


where 


T* 4 = ^T 4 


B u <Xv 


The problem, with flow properties being known across the stagnation line, is now 
reduced to the numerical evaluation of the integral 


- ‘iih 


for each spectral step in the absorption coefficient model. The optical thickness r 
given by the integral 


ri ' C( pk *) dy 


is evaluated with the aid of closed Adams quadrature formulas. Hence the integral of 
equation (C8) is determinant at each point of the shock layer. With dT* 4 /dr evaluated 
as a local constant, equation (C8) integrated by parts, and the result substituted in equa- 
tion (C5), the following expression for the divergence of the radiation flux vector at the 
kth mesh point results: 


F R,k = 4 Wk T k 4 - 2 1 Wk 2T f, 4 k - T u e 2(v) 


Mrr ?,.^ 4 - (t* 


- Vl)- 


(hi " T bk|) 
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with 



(Cll) 


where the shock layer contains M + 1 mesh points. 

A considerable saving in computing time and computer storage requirements can 
be effected if the shock layer is divided into equal increments of the optical thickness 
rather than Y. This division is carried out for each spectral step with third- order 
Lagrangian extrapolation polynomials used to obtain flow quantities at the resulting mesh 
points. The quantity Atj now becomes 


At i = 


r i.M+l 


(C12) 


For purposes of accuracy, a finer mesh is used in the calculation of F^ ^ than is used 
in the general flow-field calculation. 
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STEP MODEL ABSORPTION COEFFICIENT 

The absorption coefficient models used in this report are either identical to, or 
derived from, the model absorption coefficient reported by Olstad (ref. 8). In refer- 
ence 8, Olstad presents data, tabulated for a range of temperatures and densities for each 
step of a nine- step model of the absorption coefficient. This information, for purposes 
of reference is presented herein as table I. Included in this absorption coefficient model 
are contributions due to atomic line transitions, N + - and 0 + - electron recombination, free- 
free transitions, and molecular bond transitions. Equilibrium air properties are assumed 
to exist. 

The absorption due to line radiation in the wavelength interval 911 to 1800 A was 
determined with the aid of parameters presented by Griem (ref. 11). In order to incor- 
porate the resultant absorption coefficient information realistically into a step model, 
Olstad developed a "picket fence" model in which the line information is represented by a 
large number of narrow (in wavelength) steps of height 0 !i(p,T) superposed on a gray 
background of height q! 2 (p,T). The relative width and spacing of the steps are taken into 
account by a parameter /3(p,T). Sketch D1 illustrates Olstad’s model for line radiation 
between 911 and 1800 A. 



A detailed discussion of the model, its restrictions, and the means by which the 
governing parameters oi\, a 2 , and are determined are presented in reference 8. 

It suffices here to say that this model appears to be a reasonable means of treating the 
absorption coefficient due to line radiation where the lines may not be optically thin. 
Above 1800 A, contributions to the absorption coefficient, including line radiation, are 
treated as optically thin and are represented by a Planck average. In fact, for some of 
the larger nose radii, the line centers may be optically thick not thin. The spectral dis- 
tribution of the black body function coupled with arguments in the section "Results and 
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Discussion" (regarding the relative contributions of line centers and line wings) indicate 
that the resulting heat-transfer error will be small. 

In the vacuum ultraviolet, contributions to the absorption of radiation due to N + - 
and 0 + -electron recombination were obtained from Hahne (ref. 17). Of the twelve photo- 
ionization edges indicated by Hahne, four were used in the construction of the step model 
absorption coefficient presented in reference 8. These edges are located at wavelengths 
of 852, 911, 1020, and 1130 A. No contributions to the absorption coefficient are included 
below 400 A. 

Absorption coefficient models used in this report include a gray, a two-step, a 
three- step, and a nine- step model. The nine- step model is used directly as presented in 
table I. The other models are constructed from the nine- step model by means of Planck 
averages of various combinations of steps and are summarized in table II. 

It should be noted that all absorption models are constructed from the same basic 
spectral information and hence are consistent. This allows meaningful comparisons of 
flow profiles and heat-transfer solutions to be made between models of varying spectral 
complexity. 
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TABLE I.- NINE-STEP-MODEL ABSORPTION COEFFICIENT 
(Final numbers listed with values of pkj, pk 2 , etc., represent powers of 10; for example, 5.001+1 signifies 5.001 x 10 



BQ 

faD 




* 


s 
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b Absorption due to line steps, free-free, and free-bound interactions; 
C A11 absorption effects. 




TABLE II.- SUMMARY OF STEP- MODEL ABSORPTION COEFFICIENT 




Wavelength 

Formed from Planck 

Model 

Step 

interval, 

average of these steps 



A 

of table I 

Gray 

1 

400 to 100 000 

1 to 9 

Two step 

1 

400 to 1 130 

1 to 6 

2 

1 130 to 100 000 

7 to 9 

Three step 

1 

400 to 1 130 

1 to 6 

model A 

2 

1 130 to 1 800 

8 


3 

1 130 to 100 000 

7 and 9 

Nine step 

The nine- step model is used as it appears in table I 
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TABLE IE.- SUMMARY OF MOLE FRACTIONS OF SELECTED AIR SPECIES AT 
SHOCK WAVE AND BODY FOR V*, = 15 240 m/s, P K jp 0 = 10" 4 , and R N = 1 m 



Mole fraction 

Species 

Shock wave 

(T = 13 750 °K; p = 0.264 atm) 

Body 

(T = 8 700 °K; p = 0.275 atm) 

E" 

3.52 x 10-1 

1.14 x 10-2 

N 

2.24 x 10-1 

7.81 x 10-1 

N+ 

2.97 x 10" 1 

9.66 x lO- 3 

0 

7.00 X 10" 2 

1.87 x 10-1 

0+ 

5.31 x IQ" 2 

1.66 x lO- 3 

N2 

1.75 x 10"6 

5.13 X lO" 3 

N2+ 

2.28 x 10"6 

2.56 x 10-5 

02 

7.68 x 10“ 9 

8.33 x 10" 7 

02+ 

3.78 x 10-8 

1.10 x 10- 7 

02- 

2.35 X 10-13 

3.88 x 10-12 

NO 

2.00 x 10" 7 

7.81 X 10-5 

N0+ 

2.54 x IQ -6 

9.47 x lO" 5 
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Nondimensional coordinate 



Shock-layer coordinate, Y 


Figure 4.- Typical flow profiles. Nongray, nonadiabatic, absorbing flow; Voo= 14 000 m/s; Rjg = 10 m- p^/pg = 10~ 4 ; three-step model. 


Nondimensional coordinate 


nthalpy ratio, h/h 











L°g lf) of the nondimens ional radiative heat flux 


Total 



A 4 r3r6 rb r% 


^ NINE STEP 


imp 


A — J- 

Step 

Wavelength 


Interval, a 

1 

400 to 852 

2 

852 to 911 

3,4 

911 to 1020 

5,6 

1020 to 1130 

7,8 

1130 to 1800 

9 

1800 to 100 000 


Nose radius, R^, m 


Figure 8.- Contribution of various spectral regions to total radiative heat transfer. = 15 240 m/s ; p o 

nine-step- model absorption coefficient. 






Nondimensional enthalpy ratio, h/h 



Figure 9.- Effect on enthalpy profile of step model variations. V 00 = 16 000 m/s ; R N = 8.0 m ; p„p 0 



Nose radius, R^, m 


Figure 10.- Comparison of stagnation -point heat-transfer rates. Z = 67.06 km. 
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Nondimensional stagnation heat-transfer rate. 



Velocity, V ro> m/s 


Figure 11.- Comparison of stagnation-point heat-transfer rates. 6 a = 0.10 m. 
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Nose radius, % 


(a) Z = 42.67 km. 

Figure 12.- Variation of q R with Rf^, Voo, and Z. Calculations carried out by using three-step model. 
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Nondimens ional stagnation heat-transfer rate. 



(b) Z = 48.77 km. 
Figure 12.- Continued. 
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Nose radius , R^ 


(0 Z = 54.86 km. 
Figure 12.- Continued. 
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Nose radius, R„ 
N 


(d) Z = 60.96 km. 
Figure 12.- Continued. 
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(e) Z = 67.06 km. 
Figure 12.- Continued. 
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Nondimensional stagnation heat- transfer rate. 
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Nose radius , 


(f) Z = 73.15 km. 
Figure 12.- Concluded. 
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Figure 13.- Cooling-factor correlations. 
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Figure 13.- Continued. 
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Wall heat flux 
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Figure 13.- Concluded. 
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